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Observations  made  during  the  last  ten  years  with  the  Chandra  X-ray 
Observatory  have  shed  much  light  on  the  cooling  gas  in  the  cen¬ 
ters  of  clusters  of  galaxies  and  the  role  of  active  galactic  nucleus 
(AGN)  heating.  Cooling  of  the  hot  intracluster  medium  in  cluster 
centers  can  feed  the  supermassive  black  holes  found  in  the  nuclei  of 
the  dominant  cluster  galaxies  leading  to  AGN  outbursts  which  can 
reheat  the  gas,  suppressing  cooling  and  large  amounts  of  star  forma¬ 
tion.  AGN  heating  can  come  in  the  form  of  shocks,  buoyantly  rising 
bubbles  that  have  been  inflated  by  radio  lobes,  and  the  dissipation 
of  sound  waves. 
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Introduction 

handra  observations  of  cavities  in  cluster  centers  in  the 
last  decade  have  brought  the  issue  of  feedback  to  the  cen¬ 
ter  stage  in  astrophysics.  Feedback  from  AGN  driven  by  su¬ 
permassive  black  holes  in  the  cores  of  massive  galaxies  is  now 
seen  as  a  necessary  ingredient  to  adequately  describe  the  for¬ 
mation  and  evolution  of  galaxies  and  the  clusters  in  which  they 
sometimes  reside.  The  energy  injected  by  AGN  is  required  to 
produce  a  galaxy  luminosity  function  in  simulations  (e.g.  [1]) 
that  matches  observations.  This  feedback  may  set  the  upper 
limit  to  the  observed  masses  of  galaxies  and  contributes  to 
cluster  preheating  and  the  observed  “entropy  floors”  seen  in 
cluster  centers.  It  may  affect  scaling  relations,  such  as  the 
relationship  between  X-ray  luminosity  and  temperature  and 
potentially  affects  cluster  properties  that  are  used  for  con¬ 
straining  cosmological  models,  such  as  the  gas  mass  fraction 
(e.g.  [2]).  Data  from  the  Chandra  X-ray  Observatory  over  the 
last  ten  years  reveal  AGN  at  work  in  the  centers  of  clusters, 
inflating  bubbles  that  rise  buoyantly  through  the  intraclus¬ 
ter  medium  (ICM),  sometimes  producing  shocks  and  sound 
waves. 

Many  of  the  most  spectacular  Chandra  images  of  AGN 
feedback  are  of  interactions  in  the  centers  of  cooling  flow 
(sometimes  referred  to  as  “cool  core”)  clusters  of  galaxies. 
Under  the  assumption  of  hydrostatic  equilibrium,  the  X-ray- 
emitting  gas  in  a  cluster  is  the  most  dense  at  the  cluster  cen¬ 
ter.  The  gas  will  cool  as  it  radiates,  and  the  cooling  time  is 
inversely  proportional  to  the  density  of  the  gas.  Therefore, 
gas  will  cool  at  cluster  centers  before  it  cools  at  larger  radii. 
If  the  cooling  time  is  shorter  than  the  cluster  age,  the  cen¬ 
tral  gas  will  cool,  lose  pressure,  and  outer  gas  must  flow  in  to 
maintain  hydrostatic  equilibrium  (see  [3]  for  a  review).  Data 
from  earlier  X-ray  observatories  such  as  Einstein  [4],  ROSAT 
[5],  and  ASCA  [6]  were  used  to  give  estimates  of  the  gas  inflow 
rates  of  typically  hundreds  or  thousands  of  solar  masses  per 
year.  The  classic  “cooling  flow  problem”  has  been  that  suf¬ 
ficient  quantities  of  cool  gas  and/or  star  formation  were  not 
detected  to  match  the  predictions  from  the  X-ray  data. 

Chandra  and  XMM-Newton  observations  have  shown  that 
while  the  gas  temperature  drops  in  the  centers  of  cooling 
flow  clusters,  the  temperature  reaches  a  threshold  value  below 
which  there  is  little  cooling.  This  is  seen  in  the  temperature 


profiles  as  well  as  high  resolution  spectroscopy.  An  important 
early  result  from  XMM-Newton  high  resolution  spectra  from 
cooling  flows  was  that  the  emission  lines  from  cool  gas  were 
not  present  at  the  expected  levels,  and  the  spectra  were  well- 
fitted  by  a  cooling  flow  model  with  a  low  temperature  cutoff 
[7,  8].  This  cutoff  is  typically  one-half  to  one-third  of  the  av¬ 
erage  cluster  temperature,  and  is  usually  not  much  below  1 
keV.  The  gas  is  either  then  cooling  non-radiatively  or  being 
heated  to  higher  temperatures.  Throughout  this  paper,  we 
use  the  term  “cooling  flow”  to  indicate  clusters  with  central 
drops  in  temperature  and  large  central  density  peaks  above 
that  of  a  typical  /3-model. 

There  is  a  clear  association  with  the  presence  of  a  cool¬ 
ing  flow  in  a  cluster  and  the  central  galaxy  hosting  a  radio 
source.  Using  Einstein  data,  [9]  found  as  70%  of  cooling  flow 
clusters  had  giant  elliptical  galaxies  that  emitted  in  the  radio, 
as  compared  to  only  as  20%  of  central  galaxies  in  non-cooling 
flows.  Recent  studies,  depending  on  the  sample,  find  a  frac¬ 
tion  of  up  to  100%  for  cooling  flows  and  45%  for  non-cooling 
flows  [10]  using  Chandra  data.  This  association  is  consistent 
with  the  cooling  gas  supplying  fuel  to  the  central  black  hole, 
leading  to  the  AGN  outburst.  As  seen  in  Chandra  images,  the 
radio  jets  and  lobes  associated  with  the  AGN  significantly  af¬ 
fect  the  X-ray  gas,  inflating  “cavities”  or  “bubbles”  in  the 
cluster  centers.  Clear  correspondence  between  regions  of  ra¬ 
dio  emission  and  deficits  in  the  X-ray  are  seen  in  a  number 
of  cases  (e.g.  Hydra  A  [11,  12],  Perseus  [13],  Abell  2052,  [14], 
Abell  2597  [15,  16],  Abell  2199  [17],  Centaurus  [18],  Abell  4059 
[19,  20],  Abell  262  [21,  22],  MS0735.6+7321  [23],  M87/Virgo 
[24]).  The  X-ray  gas,  in  turn,  affects  the  radio  lobes  by  con¬ 
fining  and  distorting  them.  Energy  can  be  transported  to 
the  ICM  through  bubbles  that  rise  buoyantly  to  larger  clus¬ 
ter  radii,  as  well  as  through  weak  shocks  and  sound  waves 
generated  by  AGN  outbursts.  The  energy  transferred  to  the 
cluster  gas  can  diminish  the  cooling,  giving  an  explanation  for 
the  lack  of  copious  quantities  of  cool  gas  and  star  formation 
seen  in  the  centers  of  cooling  flow  clusters.  Composite  images 
of  A2052  and  A262  are  shown  in  Figs.  1  and  2,  respectively, 
where  blue  is  the  X-ray  emission  as  seen  with  Chandra,  red  is 
optical  emission,  and  yellow/green  (in  the  cluster  centers)  is 
1.4  GHz  radio  emission. 


Reserved  for  Publication  Footnotes 


www. pnas.org/cgi/doi/10.1073/pnas. 0709640104 


PNAS  |  Issue  Date 


Volume  |  Issue  Number  |  1-6 


Report  Documentation  Page 

Form  Approved 

OMB  No.  0704-0188 

Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 

VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 

1.  REPORT  DATE 

2Q  jq  2.  REPORT  TYPE 

3.  DATES  COVERED 

00-00-2010  to  00-00-2010 

4.  TITLE  AND  SUBTITLE 

AGN  feedback  in  clusters  of  galaxies 

5a.  CONTRACT  NUMBER 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6.  AUTHOR(S) 

5d.  PROJECT  NUMBER 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Naval  Research  Laboratory ,4555  Overlook  Avenue 

SW, Washington, DC, 20375 

8.  PERFORMING  ORGANIZATION 

REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

10.  SPONSOR/MONITOR'S  ACRONYM(S) 

11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited 

13.  SUPPLEMENTARY  NOTES 

14.  ABSTRACT 

15.  SUBJECT  TERMS 

16.  SECURITY  CLASSIFICATION  OF:  17.  LIMITATION  OF 

_ _ _  ABSTRACT 

18.  NUMBER  19a.  NAME  OF 

OF  PAGES  RESPONSIBLE  PERSON 

a.  REPORT  b.  ABSTRACT  c.  THIS  PAGE  Same  OS 

unclassified  unclassified  unclassified  Report  (SAR) 

li 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std  Z39-18 


Cavity  properties 

The  sizes  of  the  cavities  or  bubbles  created  by  the  radio  lobes 
evacuating  regions  of  the  ICM  vary  widely  from  a  few  kpc  (e.g. 
Abell  262  [21,  22])  to  hundreds  of  kpc  (e.g.  MS0735. 6+7321 
[23])  in  diameter.  The  cavities  are  often  surrounded  by  bright 
rims  or  shells,  where  the  gas  from  the  cavity  interiors  has  been 
compressed.  In  addition  to  the  visual  impression  that  the  ra¬ 
dio  lobes  are  removing  the  X-ray  gas  from  the  cavities,  we  can 
quantitatively  confirm  this  by  measuring  the  mass  of  gas  in 
the  shells  and  comparing  it  to  the  mass  of  gas  that  should  fill 
the  central  region  of  the  cluster  estimated  by  extrapolating 
the  density  profile  of  the  gas  at  outer  radii  inward.  In  Abell 
2052  [14,  25],  for  instance,  the  mass  in  the  shells  compressed 
by  the  two  radio  lobes  is  approximately  5  x  1010Mq,  consis¬ 
tent  with  the  mass  predicted  to  have  filled  the  cavities,  based 
on  the  extrapolation  of  the  density  profile. 

The  X-ray  gas  in  the  compressed  shells  surrounding  the 
radio  lobes  is  generally  found  to  be  cool,  and  does  not  show 
evidence  of  strong  shock  heating  as  had  been  predicted  with 
early  models  of  radio  source  /  ICM  interaction  (e.g.  [26]). 
Cool  gas  is  seen  surrounding  the  inner  radio  lobes  of,  for  ex¬ 
ample,  Hydra  A  [27],  Perseus  A  [28],  3C  317  in  Abell  2052 
[14,  25],  and  B2  0149+35  in  Abell  262  [21].  The  cooling  time 
within  the  shells  is  typically  longer  than  the  radio  source  age, 
suggesting  that  the  gas  in  the  shells  had  to  cool  closer  in  to  the 
cluster  center  before  being  pushed  out  by  the  radio  lobes.  In 
the  case  of  Abell  2052,  the  shell  cooling  time  is  approximately 
3  x  108  yr  while  the  age  of  the  radio  source  is  approximately 
10‘  yr.  The  bright  X-ray  shells  are  sometimes  coincident  with 
optical  line  emission,  such  as  Ha,  representing  gas  that  has 
cooled  to  temperatures  of  approximately  104  K  [14,  21,  29]. 
Therefore,  while  most  of  the  cluster  gas  cools  to  only  one- 
third  to  one-half  of  the  average  ICM  temperature,  at  least 
some  quantity  of  gas  can  cool  to  much  lower  temperatures 
and  sometimes  results  in  star  formation  [30]. 

Assuming  some  geometry  for  the  shells  surrounding  the 
radio  lobes  and  fitting  the  spectra  from  these  regions,  the 
densities  and  pressures  can  be  determined.  Alternatively,  the 
pressure  in  the  X-ray  gas  is  determined  from  undisturbed  re¬ 
gions  in  the  cluster  at  the  same  radii  as  the  cavities,  yielding 
similar  results  [31].  The  pressures  found  are  often  higher  (by 
approximately  an  order  of  magnitude)  than  the  pressures  de¬ 
rived  from  radio  observations,  assuming  equipartition,  for  the 
cavities.  The  cavities  would  be  expected  to  collapse  quickly  if 
they  were  indeed  underpressured,  and  it  is  therefore  assumed 
that  some  additional  source  of  pressure  must  be  filling  the  cav¬ 
ities.  The  origin  of  this  additional  pressure  is  still  unknown, 
but  candidates  are  low  energy,  relativistic  electrons,  a  higher 
ratio  of  ions  to  electrons  than  is  typically  assumed,  a  magnetic 
field  with  a  value  higher  than  the  equipartition  value,  or  very 
hot,  diffuse,  thermal  gas.  Hot,  diffuse,  thermal  gas  has  not 
been  detected  filling  the  radio  cavities,  but  limits  have  been 
placed  in  a  few  cases  (>  15  keV  in  Hydra  A  [27],  >11  keV  in 
Perseus  [28],  and  >  20  keV  in  Abell  2052  [25]). 

The  cavities  are  eventually  expected  to  break  up  due  to 
Rayleigh- Taylor  and  Kelvin-Helmholtz  instabilities.  Rayleigh- 
Taylor  instabilities  are  a  possible  explanation  for  the  spur  of 
gas  jutting  into  the  N  cavity  in  Abell  2052  [14,  25,  32].  The 
cavities  can  be  maintained  against  these  instabilities  by  mag¬ 
netic  draping  [33]  or  pressure  related  to  bubble  inflation  [34]. 


Buoyant  bubbles 

Since  the  density  inside  the  radio  cavities  is  much  lower  than 
that  of  the  ambient  gas,  the  cavities  should  be  buoyant,  and 
will  rise  out  to  larger  cluster  radii.  Such  outer  cavities  are 


called  “ghost  cavities”  or  “ghost  bubbles”  because  they  are 
thought  to  be  the  result  of  earlier  AGN  outbursts.  In  addition, 
if  these  cavities  are  filled  with  radio  lobes,  the  lobes  often  have 
steep  spectra,  and  radio  emission  at  the  commonly-observed 
radio  frequency  of  1.4  GHz  is  sometimes  faint  or  absent.  Radio 
observations  at  lower  frequencies  have  shown  that  the  radio 
sources  extend  into  these  outer  cavities  in  many  cases. 

Examples  of  systems  with  outer  cavities  are  Perseus  [13], 
Abell  2597  [15,  16],  Abell  262  [22],  Abell  2052  [32],  and  Hydra 
A  [12].  In  all  of  these  cases,  when  the  cavities  seemed  devoid 
of  radio  emission  at  higher  frequencies,  lower  frequency  radio 
emission  was  found  to  extend  farther  from  the  cluster  center 
to  fill  the  outer  cavities. 


Entrainment  of  cool  gas.  Cool  gas  can  be  entrained  along  with 
radio  lobes  as  they  propagate  into  the  ICM.  There  is  a  clear 
association  with  arcs  of  cool  gas  and  the  radio  jets  and  lobes  in 
the  M87/ Virgo  system  [24,  35,  36].  This  is  also  seen  in  the  cen¬ 
ter  of  Abell  133,  where  a  cool  X-ray  plume  appears  related  to 
a  detached  radio  lobe  or  relic  [37].  In  the  Perseus  cluster  [29], 
plumes  of  cool,  ss  104  K  gas  seen  in  Ha  appear  to  have  risen 
away  from  the  cluster  center,  having  been  dragged  by  one  of 
the  outer  cavities.  The  metallicity  structure  is  also  consistent 
with  the  entrainment  scenario  [38] .  Analysis  of  deep  Chandra 
observations  of  the  Hydra  A  cluster  reveals  higher  metallicity 
gas  along  the  radio  jets  and  lobes,  indicating  that  the  gas  is 
being  redistributed  from  the  center  (where  the  gas  metallicity 
is  typically  higher)  to  larger  radii  in  the  cluster  [38]. 


X-ray  cavities  as  radio  calorimeters.  Rather  than  relying  on 
radio  observations  to  estimate  the  energy  input  from  the  ra¬ 
dio  sources  where  the  ratio  of  radiative  power  to  kinetic  power 
can  vary  widely,  we  can  instead  use  observations  of  the  effects 
that  the  radio  sources  have  on  the  X-ray-emitting  ICM  to  de¬ 
termine  how  much  energy  is  being  deposited  into  the  cluster 
gas.  As  in  [11]  and  [39],  the  energy  input  is  a  combination  of 
the  internal  energy  in  the  cavity  and  the  work  done  in  inflating 
the  cavity,  or  the  enthalpy: 
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where  P  is  the  pressure  measured  in  the  X-ray  shells  (which 
are  assumed  to  be  in  pressure  equilibrium  with  the  cavities), 
and  V  is  the  volume  of  the  cavities.  The  value  of  7  is  the 
mean  adiabatic  index  of  the  gas  filling  the  cavities  and  is  4/3 
for  relativistic  gas  or  5/3  for  non-relativistic,  monatomic  gas. 
In  the  relativisitic  case,  the  total  energy  input  is  then  4 PV, 
and  in  the  non-relativistic  case  it  is  (5/2 )PV.  To  estimate  an 
average  energy  input  rate,  we  need  the  radio  source  outburst 
rate.  In  some  cases,  this  can  be  measured  directly  from  the  X- 
ray  data,  when  there  are  multiple  sets  of  bubbles.  Taking  the 
projected  distances  between  outer  cavities  and  inner  cavities 
and  calculating  the  time  to  rise  buoyantly  to  the  outer  posi¬ 
tion,  the  time  between  outbursts  can  be  estimated.  The  time 
between  outbursts  is  typically  a  few  tens  of  millions  of  years 
(e.g.  Perseus  [13];  A262  [22],  A2052  [32]).  In  cases  without 
outer  cavities,  an  approximate  value  based  on  other  systems 
is  assumed. 


Does  energy  input  from  bouyantly  rising  bubbles  offset  cool¬ 
ing?.  We  can  compare  the  energy  injection  rate  from  AGN- 
inflated  cavities  that  rise  buoyantly  in  clusters  to  the  lumi¬ 
nosity  of  cooling  gas  in  the  cluster  centers.  This  cooling  lumi¬ 
nosity  is  gvien  by 


Lcool  =  \—kT 
2  fj,mp 


[2] 
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where  M  is  the  mass  deposition  rate  into  the  cluster  center  of 
the  cooling  gas,  T  is  the  upper  temperature  from  which  the 
gas  is  cooling,  and  /x  is  the  mean  molecular  weight  with  fj,mp 
being  the  average  mass  per  particle  in  the  gas.  In  many  cases, 
the  energy  input  rate  from  the  buoyant  bubbles  is  enough  to 
offset  the  cooling  gas  (e.g.  Abell  2052  [14,  25,  32]  with  a  rate 
of  «  3  x  1043  erg  s-1  and  Hydra  A  [11,  27,  40]  with  a  rate 
of  «  3  x  1044  erg  s_1).  In  some  cases,  the  bubbles  did  not, 
at  least  initially,  appear  to  supply  enough  energy  to  the  ICM 
(e.g.  Abell  262  [21],  where  the  energy  input  rate  is  approxi¬ 
mately  an  order  of  magnitude  too  low). 

In  addition  to  detailed  studies  of  individual  objects,  stud¬ 
ies  of  samples  of  cooling  flow  clusters  with  central  AGN  have 
been  undertaken  [41,  42,  43,  44,  45].  Thirty-three  systems 
were  included  in  a  sample  in  [45],  and  it  was  found  that,  on 
average,  4 PV  of  enthalpy  per  cavity  supplies  enough  energy 
to  substantially  suppress  cooling  in  clusters. 


Shocks 

Simulations  of  AGN  outbursts  in  cluster  centers  show  the  cre¬ 
ation  of  shocks  as  well  as  cavities  related  to  the  radio  lobes 
(e.g.  [46]).  These  “cocoon”  shocks  can  heat  the  cluster  gas 
not  only  in  the  direction  parallel  to  the  propagation  of  the  ra¬ 
dio  lobes,  but  perpendicular  to  their  propagation  axis,  as  well. 
Therefore,  energy  can  be  distributed  throughout  the  cluster 
center  if  shocks  are  present. 

Until  recently,  shocks  were  not  seen  associated  with  AGN 
radio  lobes  in  cluster  cooling  flows,  but  in  the  last  few  years, 
several  have  been  detected.  A  spectacular  example  is  the  clus¬ 
ter  MS0735. 6+7321  [23]  which  hosts  extremely  large  cavities 
with  diameters  of  approximately  200  kpc.  The  total  energy 
injection  required  to  inflate  the  cavities  and  produce  the  ob¬ 
served  shocks  is  6  x  1061  erg,  making  this  the  most  powerful 
radio  outburst  known.  The  shock  energy  input  rate  is  given 
by 
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where  P  is  the  pre-shock  pressure,  7  =  5/3,  and  2tt/uj  is  the 
time  interval  between  shocks  [31].  A  deep  Chandra  image  of 
M87/Virgo  revealed  an  approximately  spherical  shock  in  the 
cluster  center  [24].  The  shock  is  most  clearly  seen  in  a  hard- 
band  image,  and  is  fairly  weak  with  a  Mach  number  of  1.2. 
A  Mach  1.65  shock  was  found  in  the  Hercules  cluster  [47], 
with  the  total  energy  deposited  3  x  1061  erg.  There  is  a  weak 
shock  (Mach  1.3)  outside  of  one  of  the  bubble  rims  in  the 
center  of  the  Perseus  cluster  [48].  As  is  common,  the  density 
jump  measured  is  consistent  with  a  shock  model,  but  a  rise 
in  temperature  is  not  measured.  Two  shocks  are  seen  sur¬ 
rounding  the  center  of  Abell  2052  [32],  both  with  Mach  1.2. 
The  temperatures  across  the  shocks  are  consistent,  within  the 
errors,  to  the  jumps  expected  in  the  shock  models,  but  the 
best-fitting  temperature  values  are  constant  across  the  fronts. 
If  an  increase  in  temperature  is  not  measured  behind  a  shock, 
[49]  has  suggested  that  shocks  may  be  isothermal  if  conduc¬ 
tion  is  efficient  in  these  regions.  However,  the  temperature 
rises  are  very  difficult  to  detect  given  the  narrow  width  of  the 
shock  fronts,  adiabatic  expansion,  and  the  difficulty  in  mod¬ 
eling  the  large  amount  of  cluster  gas  at  larger  radii  projected 
onto  these  regions  when  fitting  the  X-ray  spectra.  Based  on 
further  analysis  of  the  Perseus  system,  [48]  conclude  that  mix¬ 
ing  of  post-shock  gas  with  cool  gas  associated  with  optical-line 
emission  may  result  in  the  lack  of  temperature  rise  seen  asso¬ 
ciated  with  the  shock. 


Sound  waves 

In  addition  to  weak  shocks  in  the  centers  of  cooling  flow  clus¬ 
ters,  sound  or  pressure  waves  have  also  been  observed.  These 
features  are  most  prominent  in  the  Perseus  cluster  [49,  50]. 
They  appear  as  5  —  10%  increases  in  pressure  in  ripple-like 
patterns  surrounding  cluster  centers.  Changes  in  tempera¬ 
ture  are  not  seen  in  these  regions,  and  it  is  likely  that  their 
energy  is  dissipated  by  viscosity.  There  may  be  one  or  more 
set  of  ripples  corresponding  to  each  AGN  outburst  depending 
on  variations  in  energy  output  that  can  occur  during  a  single 
outburst.  Assuming  they  are  associated  with  single  outbursts, 
the  ripple  separations  in  Perseus  imply  a  repetition  rate  of  10' 
yr  [50].  As  shocks  propagate  outward  from  cluster  centers, 
they  weaken  and  may  appear  as  sound  waves.  In  addition  to 
the  Perseus  cluster,  sound  wave  ripples  are  seen  in  the  Cen- 
taurus  cluster  [51]  and  M87/Virgo  [24].  Concentric  ripple-like 
features  are  also  seen  surrounding  the  center  of  Abell  2052, 
but  current  analysis  shows  that  these  are  consistent  with  weak 
shocks  [32]. 

A  key  difference  between  shock  and  sound  wave  heating 
is  that  sound  wave  dissipation  depends  on  transport  coeffi¬ 
cients  and  shock  dissipation  does  not  [31].  The  sound  wave 
dissipation  rate  is  sensitive  to  temperature  and  would  be  less 
effective  in  lower  temperature  systems.  This  may  be  offset  by 
a  higher  frequency  of  outbursts  in  cooler  clusters  [31].  One  of 
the  important  implications  of  the  sound  waves  is  that  energy 
can  be  deposited  over  4-7r  steradians,  and  not  only  along  the 
narrow  angle  of  the  radio  jets. 


Case  studies:  Abell  2052  and  Abell  262 

Abell  2052  is  one  of  the  few  clusters  that  showed  evidence 
for  radio  bubbles  using  data  from  an  X-ray  telescope  prior  to 
the  Chandra  /  XMM-Newton  era  ( ROSAT  [52]).  It  was  first 
observed  with  Chandra  for  37  ksec  in  Cycle  1  in  2000,  and 
subsequently  observed  for  129  ksec  in  Cycle  6  in  2006  and  497 
ksec  in  Cycle  10  in  2009.  A  three-color  image  from  the  total 
657  ksec  cleaned  dataset  is  shown  in  Fig.  3,  where  the  col¬ 
ors  represent  different  energy  ranges  (red  0.3  —  1.0  keV,  green 
1.0  —  2.0  keV,  and  blue  2.0—  10.0  keV).  Clear  bubbles  are  seen 
to  the  north  and  south  of  the  cluster  center,  and  a  filament  ex¬ 
tends  into  the  northern  bubble.  Surrounding  the  bubbles  are 
bright  shells  of  emission,  and  outside  of  these  shells  a  jump  is 
the  surface  brightness  is  visible  extending  around  the  cluster 
center.  A  second  jump  is  seen  outside  this  one,  and  it  appears 
sharper  to  the  NE  than  to  the  SW.  It  is  possible  that  this 
feature  is  a  cold  front  rather  than  a  shock. 

An  unsharp-masked  image  for  the  total  657  ksec  dataset 
is  shown  in  Fig.  4  in  grayscale  with  1.4  GHz  radio  emission 
from  the  VLA  FIRST  survey  [53]  superposed  in  green  con¬ 
tours.  Concentric  ripple-like  features  are  seen  surrounding  the 
cluster  center,  and  these  are  modeled  as  shocks  in  [32]  based 
on  the  earlier  163  ksec  dataset.  These  features  are  at  31  and 
46  kpc  from  the  AGN  and  the  density  jumps  associated  with 
them  are  both  consistent  with  shocks  with  Mach  number  1.2. 
The  separation  of  the  potential  shock  features  is  consistent 
with  a  radio  source  outburst  rate  of  2  x  10'  yr.  This  same 
cycle  time  is  found  when  considering  the  locations  of  outer 
cavities  in  the  system.  A  small  cavity  appears  bounded  by  a 
narrow  filament  to  the  NW  of  the  N  bubble,  and  the  S  cavity 
appears  split  into  two  cavities,  with  the  outer  one  to  the  SE. 
Assuming  the  outer  cavities  rose  buoyantly  to  their  current 
positions,  we  find  a  cycle  time  of  2  —  4  x  107  yr  depending  on 
whether  they  rose  at  0.5  times  or  1  time  the  sound  speed. 

In  [32],  we  found  a  mass  deposition  rate  of  M  =  55±4Mq 
yr-1  and  a  cooling  luminosity  of  Lcooi  =  5.4x  1043  erg  s_1.  We 
find  the  rate  of  energy  input  from  the  buoyantly  rising  bubbles 
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assuming  they  are  filled  with  relativistic  plasma  (7  =  4/3)  to 
be  3  —  6  x  1043  erg  s~4  depending  on  the  buoyancy  rise  time. 
The  input  from  shock  heating  is  only  1  x  1043  erg  s_1,  how¬ 
ever,  the  combination  of  heating  from  shocks  and  buoyantly 
rising  bubbles  can  offset  the  cooling  in  this  system. 

Abell  262  was  first  observed  by  Chandra  in  Cycle  2  in  2001 
for  30  ksec.  It  was  subsequently  observed  in  Cycle  8  in  2006 
for  112  ksec.  A  mass  deposition  rate  of  M  =  19 /jMg  yr~4 
and  a  cooling  luminosity  of  Lcooi  =  1.3  x  1043  erg  s_1  were 
found  in  [21].  Based  on  the  one  clear  bubble  seen  in  the  earlier 
dataset,  it  was  found  that  the  radio  source  energy  input  fell 
more  than  an  order  of  magnitude  short  of  offsetting  the  cool¬ 
ing  [21],  and  it  was  concluded  that  the  AGN  may  have  had  a 
more  powerful  outburst  in  the  past  so  that  cooling  could  be 
offset  on  average.  In  this  study,  a  repetition  rate  of  1  x  108  yr 
was  assumed  for  the  AGN  outbursts. 

Analyzing  the  longer  dataset,  however,  we  found  in  [22] 
that  the  radio  source  can  offset  the  cooling.  This  analysis  was 
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based  on  both  a  deep  Chandra  observation  and  low-frequency 
radio  observations.  Several  more  bubbles  were  found,  and  a 
tunnel  made  up  of  multiple  bubbles  rising  with  short  time  in¬ 
tervals  was  found  to  the  W  of  the  cluster  center.  Using  the 
multiple  sets  of  bubbles,  the  radio  source  comes  within  a  fac¬ 
tor  of  two  of  quenching  the  cooling.  Chandra  X-ray  residual 
images  superposed  with  radio  emission  at  1440  and  610  MHz 
are  shown  in  Fig.  5.  It  may  well  be  that  further  analysis 
of  deeper  multi-frequency  observations  of  objects  in  samples 
such  as  [45]  that  seem  to  have  insufficient  energy  input  rates 
would  reveal  that  the  AGN  are  indeed  powerful  enough  to 
offset  cooling. 
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Fig.  1.  Composite  X-ray  (blue,  Chandra  657  ksec  image),  optical  (red,  Digitized  Sky 
Survey),  radio  (green/yellow,  1.4  GHz  [53])  5;.6  X  4/. 7  (230  X  190  kpc)  image  of  Abell 
2052. 


Fig.  2.  Composite  X-ray  (blue,  Chandra  139  ksec  image),  optical  (red,  Digitized  Sky 
Survey),  radio  (green/yellow,  1.4  GHz  [22])  5/.6  X  4/.7  (110  X  90  kpc)  image  of  Abell  262. 


Fig.  3.  Three-color  Chandra  657  ksec  4/. 3  X  3;.7  (180  X  150  kpc)  image  of  Abell  2052 
(red  =  0.3-1  keV,  green  =  1-2  keV,  blue  =  2-10  keV). 


Fig.  4.  Unsharp-masked  657  ksec  Chandra  2r.8  X  2'  A  (115  X  95  kpc)  image  of  Abell 
2052  with  1.4  GHz  radio  contours  from  the  VLA  FIRST  survey  superposed.  Multiple  sets  of 
bubbles  are  seen,  as  well  as  ripple-like  jumps  in  surface  brightness  that  are  consistent  with  weak 
shocks. 


Fig.  5.  Residual  Chandra  139  ksec  3/.7  X  l/.5  (70  X  30  kpc)  images  (green)  of  Abell 
262  with  radio  emission  superposed  in  pink  (VLA  1.4  GHz  [top]  and  GMRT  610  MHz  [bottom], 
[22]). 
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